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Abstract 

Species expand their geographical ranges following an environmental change, long range 
dispersal, or a new adaptation. Range expansions not only alter the affected ecosystems, 
but also have a significant impact on the evolution of the expanding species. Although the 
dynamics of deleterious, neutral, and beneficial alleles have been extensively studied, the fate 
of alleles under frequency-dependent selection remains largely unexplored. The dynamics of 
cooperative alleles are particularly interesting because selection can be both frequency and 
density dependent. We find that the fate of cooperation critically depends on the coupling 
between population and evolutionary dynamics at the front of the expansion. In particular, the 
population at the front can be substantially more cooperative than the population in the bulk. 
Under certain conditions, the entire front can be taken over by cooperators, which leads to the 
splitting of the cooperative population from defectors and a substantial increase in the total 
number of cooperators in the population. Unlike previously reported surfing events for neutral 
and deleterious mutations, wave splitting can be both stochastic and deterministic. Our results 
not only provide an explanation for the maintenance of cooperation, but also elucidates the 
effect of eco-evolutionary feedback on the maintenance of genetic (or species) diversity during 
range expansions. Quite surprisingly, we find that defectors can spread much faster in mixed 
waves with cooperators than they would be able to on their own or by invading cooperators. 
This enhanced dispersal in mixed waves could counterbalance the loss of genetic diversity due 
to the founder effect for mutations under frequency-dependent selection. Although we focus on 
cooperator-defector interactions, our analysis, especially of wave splitting, could be relevant for 
other systems described by reaction-diffusion equations. 

Introduction 

Cooperation between organisms has always interested evolutionary biologists [JJ EJ |3]. On the one 
hand, cooperative interactions are widespread in living systems. Microbes cooperate to digest food, 
scavenge for scarce resources, or build a protective biofilm [3]; animals cooperate to hunt or avoid 
predation [5] ; and individual cells cooperate to enable multicellular life [6j [7] . On the other hand, 
evolutionary theory has struggled to explain the mere existence of cooperation [U [2} [3] . Although 
cooperation is beneficial to the species, it is susceptible to invasion by defectors, individuals who 
reap the benefits of cooperation without paying the costs. Defectors, having a higher relative fitness, 
are expected to take over the population leading to the demise of cooperation. The breakdown 
of cooperation often has devastating consequences such as the tragedy of the commons [8] or can- 
cer P [10]. The break down of cooperation can also be desirable, e.g., when it destroys biofilms 
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protecting pathogens from antibiotics and the immune system [4]. Understanding the evolution 
and maintenance of cooperation is therefore an important problem in evolution, ecology, economics, 
medicine, and other branches of science. 

Several mechanisms have been proposed that can stabilize cooperation against defection jH El [11] , 
including direct reciprocity [T2], [13] , group and kin selection [Ml [151 [16] , an d spatial structure [T7 [ [T8 | 
[19]. Most of the previous studies have focused exclusively on the changes in the relative frequencies 
of cooperators and defectors and neglected possible changes in the population size. The naive 
expectation, however, is that a reduced level of cooperation should lead to a lower average fitness 
of the population and, therefore, to a lower population size. Indeed, this effect of evolutionary on 
ecological dynamics has been observed in several experimental populations [201 ED [22] as well as in 
ecological public goods games [23l EH [25]. In the latter studies, the authors also found that, under 
certain conditions, cooperators are favored by natural selection at low population densities while 
defectors are favored at high population densities. This dependence of evolutionary dynamics on 
population density suggests that the low-density edges of population ranges might be conducive to 
the evolution of cooperation. The edge of an expanding population is of particular interest because 
the new territories might be colonized mostly by cooperators. 

Range expansions and range shifts are common in nature [2S1 ED ESI ESI GEIl ED E21 E3]- Ex- 
amples include the recolonizations of temperate latitudes between glaciations, the invasion of North 
American forests by the Asian long-horned beetle (Anoplophora glabripennis) , and the spread of the 
western corn rootworm (Diabrotica virgifera) in the Midwestern United States. Apart from ecologi- 
cal and sometimes economic impact, range expansions bring about significant changes in the genetic 
diversity and the evolutionary dynamics of the expanding species [SI ESI EHl ED EH1 EHl HO]. In 
particular, genetic diversity is typically lost because of the founder effect, and mutations appearing 
close to the expansion edge are very likely to reach high frequencies in the population even if they 
are neutral or slightly deleterious [311 El E21 EHl S3]. Although the fates of neutral, deleterious, and 
beneficial mutations arising at the edge of a range expansion have been extensively studied, the fate 
of mutations subject to frequency-dependent selection, e.g. encoding cooperative traits, has received 
little attention. 

Here, we study how the interplay between ecological and evolutionary processes affects the evo- 
lution of cooperation during a range expansion. To this end, we formulate a model that combines 
the effect of genetic composition on population growth and the effect of population density on the 
fitnesses of different genotypes. Our works complements two previous studies of the evolutionary 
dynamics of cooperation in spatial populations [251 HI]- I n [25], the intricate spatial patters of 
coexisting cooperators and defectors were reported for ecological public goods games. In [Hj, the 
authors found that cooperation can persist in a spatial prisoner's dilemma game because of the cyclic 
turnover of cooperators, defectors, and extinctions. Similar to our work, their results rely on the 
fact that cooperators grow faster than defectors when they are spatially isolated from each other. 

We find that cooperators are favored at the edge of an expanding population and, under certain 
conditions, cooperators can outrun defectors and spread into unoccupied territories, leaving defectors 
behind. This mechanism of maintaining cooperation might play an important role in populations that 
experience frequent disturbances (e.g. forest fires) followed by range expansions. Our results suggest 
that, in such populations, cooperators may experience selection for traits that have little effect in 
non-expanding populations. More generally, we report the splitting of a mixed two-species (or two- 
allele) traveling wave into a wave of one species followed by an invasion wave of the other species. 
This wave-splitting phenomenon might also be relevant to other processes described by traveling 
react ion- diffusion waves, which are widely used in biology and chemical kinetics [45 1 f36l |4*6] . 

When splitting does not occur, defectors slow down range expansions. The rate of spreading 
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is primarily controlled by the evolutionary dynamics at the front. Interestingly, when evolutionary 
dynamics are much slower than ecological dynamics, we find that a mixed wave of cooperators and 
defectors can experience long periods of acceleration. We also find that defectors can spread faster 
in mixed waves than it would be able in isolation or by invading cooperators. This finding could be 
important for conservation efforts to help ecosystems shift in response to a rapid climate change or 
habitat deterioration. 



To understand the fate of cooperation during range expansions, we formulate a spatial model of a 
mixed population of cooperators and defectors colonizing new territories. The theoretical studies 
of range expansions were pioneered in [IT] HS] . and a good introduction to this topic can be found 
in [15]. We first discuss range expansions of pure cooperators and then consider defectors invading 
cooperators. Although these two spreading phenomena have a similar mathematical description, 
they have rarely been studied together, because one is an ecological, and the other is an evolutionary 
process. Our model allows for changes in both population density and allele frequencies and explicitly 
includes the coupling between these two variables. For simplicity, only one-dimensional expansion 
waves are considered, which should be a good approximation when number fluctuations and the 
curvature of the wave front can be neglected. 

Ecological dynamics 

We begin by considering populations of only cooperators. Range expansions are driven by migration 
from colonized to new territories and by population growth. When migration (or dispersal) is short- 
range and isotropic, it can be approximated by a diffusion term leading to the following reaction- 
diffusion model of range expansions 



where c(t, x) is the population density at time t and spatial location x, D is the effective diffusion 
constant, and G c (c) is the per capita growth rate. 

Since any habitat has a limited carrying capacity, the per capita growth rate, G c (c), must decline 
and become negative at high population densities. Populations with cooperative growth may also 
experience reduced or negative growth rates at low population densities because, for small c, the 
probability of forming cooperating groups or the size of these groups is too small. Such nonmonotonic 
dependence of the per capita growth rate on population density is called an Allee effect and has been 
observed in different species ranging from budding yeast to desert bighorn sheep [HH [501 ETJ [52j [531 
I5H [55] . The most common model of an Allee effect assumes the following growth rate [HI [50] 



where K is the carrying capacity, c* is the Allee threshold, i.e. the minimal density required for 
populations to grow, and g c K 2 sets the overall magnitude of the per capita growth rate. One typically 
distinguishes a strong Allee effect when c* > and a weak Allee effect when — K < c* < 0. Allee 
effect is absent when c* < — K because G c (c) monotonically decreases for c > 0. Thus, equation (j2J) 
is sufficiently flexible to describe populations with and without an Allee effect. 

Even if the growth rate is negative at small densities, populations can spread into unoccupied 
territories. At t — 0, we assume that the habitat is colonized for all x < 0, but it is empty for x > 0. 



Models 




(1) 




(2) 
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After an initial transitory period, expansion waves typically move at a constant velocity and the 
density profile does not change in the comoving reference frame, i.e. the reference frame moving along 
with the expansion wave; see figure [T^b. The expansion velocity is known exactly [Io" | US ] W7 \ \5E \ I5T] 
and is given by 



2& (K - 2c*) , c* > -K/2 
2^Dg c K\c*\, c* < -K/2. 



(3) 



Coordinates in the comoving reference frame (r, () are the defined in term of (t, x) as r = t and ( 
x — v c t. For c* > —K/2, the shape of the wave profile is known exactly [33], EHJ [57] : 



K 

c(C) = (4) 

1 + g V 2D ^ 



For c* < —K/2, the front has a qualitatively similar shape with c(£) declining from 1 to as ( 

j ZZ / gcK\c*\ j- 

goes from — oo to +oo. The width of this transition scales as K ^ , and c(Q ~ e~V s s f or 
large C 03 S7J M\ • 

Evolutionary dynamics 

Similar to organisms, alleles encoding cooperator or defector phenotypes can spread in populations. 
Under the assumption of short-range and isotropic migration, the spreading of genetic changes can 
also be modeled by a react ion- diffusion equation: 

df d 2 f 

provided the alleles do not affect how organisms migrate and disperse. Here, f(t,x) G (0, 1) is the 
frequency (fraction) of defectors, and Gf(f) is the relative growth rate of defectors describing the 
force of frequency- dependent natural selection. Note that the diffusion constant in equation (J5j) is 
the same as in equation (CQ) because both genetic and population spreading are due to the same 
migration process. 

The following model of frequency- dependent selection is most commonly used because of its 
simplicity and because it appears naturally as a weak- select ion limit of evolutionary game theory [21 

E7JEB] 

G f (f)f = 9ff(l-f)(r-f), (6) 

where gj > is the strength of selection and /* is the preferred or equilibrium frequency of de- 
fectors. In spatially homogeneous populations, cooperators and defectors coexist at a stable fixed 
point f — f*, provided /* G (0, 1). The coexistence of cooperators and defectors has been observed 
experimentally; see, e.g., When /* > 1, defectors outcompete cooperators, while cooperators 
prevail when /* < 0. Negative gf and /* G (0, 1) describe a bistable behavior, e.g. due to a chemical 
warfare, and is not considered here. In game theory, these four scenarios are known as the snowdrift, 
prisoner's dilemma, harmony, and coordination games [58] . 

To understand the maintenance of cooperation, we need to know how defectors invade coopera- 
tors; see figure [Tbd. The velocity of this invasion can be calculated using the results of [IS] and is 
given by 
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Vi = 2y / J DG / (0) = 2 V AD^7 7 . (7) 

In the comoving reference frame, the frequency of defectors changes from to /* (or from to 1, 
if /* > 1) as ( goes from — oo to +oo. The characteristic width of this transition scales as \J D /(gff*). 

Coupling between ecological and evolutionary dynamics 

In the preceding discussion of population and genetic waves, we avoided the coupling between ecology 
and evolution either by assuming a constant genetic composition (no defectors) or neglecting the 
changes in population properties (e.g. carrying capacity) due to defector invasion. More general 
situations require a combined model with G c and Gf depending on both c and /. Here, we consider 
a natural extension of equations ([SJ and (jHJ): 



G c (c,f)c = g c (f)c[K(f)-c][c-c*(f)} 

Gf(cJ)f = g f (c)f(l-f)[r(c)-f], [ } 

where the parameters of population dynamics depend on the genetic composition, and the parameters 
of evolutionary dynamics depend on population density. Two comments on the functional form of 
our model are in order. First, equation ((§]) allows for the most general dynamics that reduces to the 
classic models of frequency- dependent selection and cooperative growth with an Allee effect. Second, 
G c (c,f) and Gf(c,f) of an arbitrary functional form can be recast in the form of equation (jHJ) by 
allowing g c to depend on c and gf on f. We expect these dependences to be small in many situations 
because the other terms in equation (JH]) describe the most important aspects of the population 
dynamics. The analysis presented below, however, does not depend on the assumption that g c is 
only a function of / and gf in only a function of c, provided g c (c, f) > 0, <?/(c, /) > 0, and Gf is a 
decreasing function of / for / e (0, /*). In the following, we will assume that G c and Gf satisfy these 
conditions and will illustrate our results in the context of a simpler model given by equation (jSJ). 
The advantage of this approach is that many calculations can be carried out explicitly in the simpler 
model thus allowing us to provide an intuitive interpretation of the results. 

The behavior of g c (f) and <7/(c) has not been previously investigated; therefore, we typically 
assume that these two functions are constants. The dependencies of other model parameters in 
equation (jSJ) on c or / are known qualitatively from previous experimental and theoretical studies. 

Several experiments have confirmed the naive expectation that K(f) is a decreasing function [2H 
|2"U| 122] . Interestingly, the opposite effect of defectors has also been observed in recent experiments 
with budding yeast, where it was established that mixed populations of defectors and cooperators 
have a higher carrying capacity than pure populations of cooperators [60] . For simplicity, we assume 
that K is a constant in most numerical solutions. 

The Allee threshold c* is expected to increase with /, and this was indeed observed both in 
models [23] and recent experiments [HI]. In numerical solutions, we use c*(f) = c*(0)/(l — /), which 
is equivalent to requiring a certain density of cooperators to produce a sufficient amount of public 
goods for population growth. 

As we show below, the dependence of the preferred frequency /* on population density is particu- 
larly important for the fate of cooperation during range expansions. Modeling of public goods games 
revealed that, under certain conditions, /* is a increasing function of c [231 [21], which has recently 
been confirmed by experiments with budding yeast [61]. For simplicity, in numerical solutions we 
assume that /* changes from a high to a low value at a critical population density c; in other words, 
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we use /*(c) = /*(0) + [/*(1) - /*(0)] i?(c - c), where #(•) is the Heaviside step function, which 
equals one for positive arguments and zero for negative arguments. Note that the value of /(0) can 
be negative when cooperators outcompete defectors at low population densities. 

A spatial model with G c and Gf defined in equation (jSJ) cannot be obtained by simply combining 
equations (JTJ and (jSJ) because changes in population density and defector frequency due to migration 
are not independent. Instead, migration terms have to be added to the dynamics of cooperator and 
defector densities, defined cLS Cq — C (1 — /) and Cd = cf. Upon using equation (jSJ) to calculate the 
growth rates for c c and Cd and adding the diffusion terms, we obtain that the spatial dynamics of 
mixed populations can be described by the following equation 



dc d 2 c 

-dF = D ^ + GcCc ~ GfCd 



dcd 

at 



D 



9% 
dx 2 



G c c d + G f c d . 



(9) 



Here, the terms with G c come from the ecological dynamics controlling the population density, and 
the terms with Gf come from the evolutionary dynamics controlling the relative abundances of 
cooperators and defectors. Note that we used the same diffusion constants for both cooperators 
and defectors. This is justified as long as the mutations causing defector phenotype do not affect 
dispersal. 

A range expansion of a mixed population is shown in figure [Tfef. Similar to the pure cooperator 
and invasion waves, mixed waves can spread with a constant velocity v m after an initial transient. 
Both cooperator and defector density profiles reach a constant shape in the comoving reference 
frame, but these shapes are different. In particular, the relative frequency of cooperators is much 
higher at the front than in the population bulk. This behavior is expected because /* is a decreasing 
function of c, and population density declines at the front. 

It is convenient to separate evolutionary and ecological dynamics, so we recast equation (jHJ) in 
terms of population density and defector frequency: 



Be f) 2 c 

^ = D^ + G c (cJ)c 

Of d 2 f 91 df (10) 

dt dx 2 dx dx ~^ * ' ^ ^ 

Note that the additional advection term 2D 9 || appears because changes in f(t,x) due to mi- 
gration depend on c(t,x). In particular, migration increases the relative frequency of organisms in 
low-density regions if these organisms are more abundant in nearby high-density regions because 
regions with higher density send out more migrants. This effect of density gradients is particularly 
important at the wave front, where c(t, x) is changing rapidly, and there is a constant imbalance 
between migrants coming from the high-density colonized regions and migrants coming from low- 
density uncolonized regions. A more general discussion of the advection terms in populations with 
many different species (or alleles) can be found in |36j . 



Numerical solutions and simulations 

Numerical solutions of equation ()9]) can be easily obtained by standard methods [62]. Equations (JT]) 
and (jSJ) did not require a separate solver because they are special cases of equation (jUJ). We used an ex- 
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population wave 



genetic wave 



mixed wave 




Figure 1: Population, genetic, and mixed traveling waves, (a) A population of cooperators (blue) 
expands into an empty territory (black) with a constant velocity v c . (b) shows the density profile 
of the population wave in (a), (c) Defectors invade cooperators establishing a mixed population. 
The color encodes defector frequency with blue corresponding to pure cooperators and red to pure 
defectors. Mixed populations then have a magenta color. The invasion wave moves at a constant 
velocity v\. (d) shows the defector density profile of a genetic wave in (c). (e) A mixed population 
of defectors and cooperators expands into empty territory with a constant velocity v m . Initially, 
cooperators and defectors are homogeneously mixed in equal proportion. As the population expands, 
the relative frequency of cooperators increases at the wave front, but stays constant in the interior 
of the population, (f) The density profiles of defectors and cooperators show the enhancement of 
cooperation at the expansion front and a small lag between defectors and cooperators. For all panels 
in this figure, we used K = 1, c*(f) = 0.2/(1 - /), f*(c) = 5tf(c - 0.5) - 4.5, g c = 2, g f = 0.25, 
D = 0.5, and numerically solved equation ([9]). The wave profiles are plotted for t = 94. 
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plicit forward-time centered-space (FTCS) finite-difference method (4-point stencil) [62] • Spatial dis- 
cretization step was 0.1, and temporal discretization step was min{10~ 3 /-D, 10~ 2 / max{gj max{|/*|, 1}, g c K 2 }}. 
This level of discretization was sufficient to ensure that numerical wave velocities do not differ from 
the analytical results in equations fl3]) and ([7]) by more than a percent. For cooperator and mixed 
waves, the initial population typically occupied the left 10% of the habitat. For genetic waves, the 
whole habitat was occupied by cooperators, while defectors were initially present only the left 10% of 
the habitat. At t = 0, populations were always in local equilibrium, i.e. with c = K(f) and / = f*(c). 
We used no-flux boundary conditions and computed solutions up to the time when the wave had 
spread into 90% of the habitat. The habitat size L, i.e. the distance between the left and right 
boundaries, was typically equal to 100. 

We also performed individual-based simulations to demonstrate the possibility of stochastic 
splitting. The simulations were done on a one-dimensional lattice of sites each with a carrying 
capacity K. Every time step consisted of one possible migration event and one reproduction 
or death event at every site. During a migration event, a randomly chosen organism migrated 
with probability m to one of the two neighboring sites. Migration was isotropic, and we im- 
posed no- flux boundary conditions. During a reproduction or death event, the population at a 
site could increase by one, decrease by one, or remain unchanged. The probability of birth was 
given by g c c 2 (K + c*)/K 3 /(l + g c c 2 (K + c*)/K 3 + g c (c 3 + cKc*)/K 3 ), and the probability of death 
was given by g c (c 3 + cKc )/K 3 /(l + g c c 2 (K + c*)/K 3 + g c (c 3 + cKc*)/K 3 ). Death events were 
equally likely to eliminate cooperators and defectors, but birth events created either a new cooper- 
ator or a new defector with probabilities determined by the interaction matrix A which was chosen 
to mimic equation (jSJ). Given that a birth event occurred, a cooperator was born with probabil- 
ity f(A cc (l - /) + A cd f)/(A CC (1 - f) 2 + A cd (l - f)f + A dc f(l -f)+ A dd / 2 ), where / is the frequency 
of defectors, and we used letters c and d as indices. Otherwise, a defector was born. K time steps 
constituted a generation. 

Results 

In the previous section, we formulated a reaction-diffusion model that includes the coupling between 
ecological and evolutionary dynamics of cooperators and defectors. This model, equation (fTUj) . and 
its special cases, equations (pQ) and (jSJ), describe range expansions of cooperators with velocity v c , 
invasion of cooperators by defectors with velocity V{, and spreading of mixed waves of defectors and 
cooperators with velocity v m . For mixed waves, we find that the frequency of cooperators is higher 
at the wave front because cooperators are favored at low population densities; see figure [Tfef. In this 
section, we show that, under certain conditions, cooperators can take over the entire wave front and 
split from defectors by colonizing new territories faster than defectors can invade from behind. We 
also investigate how the speed of mixed waves depends on the parameters of the model and show 
that, when evolutionary dynamics are much slower than ecological dynamics (<?/ <C g c K 2 ), mixed 
waves can experience long periods of acceleration. 

The behavior of mixed waves depends on the ratio of cooperator and invasion velocities. When v ; > 
v c , defectors invade faster than cooperators can spread into new territories; therefore, any initial 
condition leads to a mixed population of defectors and cooperators. This mixed population will ex- 
pand into empty territories with both defectors and cooperators spreading at the same velocity v m . 
When v c > v- u two qualitatively different outcomes are possible. Either cooperators and defectors 
can spread together in a mixed wave with velocity v m (figure [2^bc), or a mixed wave can split into an 
ecological expansion wave and an evolutionary invasion wave (figure EUef) . In the latter scenario, the 
population expansion with velocity v c is driven solely by cooperators followed by a slower invasion 
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by defectors with velocity v- x . 

The analysis of equation (TlOl) is complicated both by the coupling between the two differential 
equations and by the unknown dependences of the parameters on population density and defector 
frequency. To reduce this complexity, we first consider the effect of evolution on ecology and the 
effect of ecology on evolution separately and then analyze the general case. 



Effects of evolutionary on ecological dynamics 

By neglecting the effect of ecology on evolution and setting <?/(c) = gj and /*(c) = /*, we can imme- 
diately see that f(t,x) = f* is a stationary solution of equation ( flOl) . With this solution for f(t,x), 
the dynamics of c(t, x) become identical to that in equation ([1]) with the parameters c*, K, and g c 
evaluated at / = /*. Therefore, the velocity and density profile of the mixed wave can be immedi- 
ately obtained from equations flSJ) and (jl]), e.g. the mixed velocity v m is given by 

= { y ^F [K(n-2c*(n] , c *(f*) > -k/2 

1 2^/Dg c (f*)K(f*)[-c*(f*)], c*(/*) < -K/2. 

Defectors are expected to decrease the rate of spatial expansions [63J. Consistent with this expecta- 
tion, equation (fTTj) predicts that Vf* < v c , provided that all or some of the growth parameters (g c , K, 
and c*) substantially decrease with /. Note that, because /* is an attracting fixed point everywhere 
in space, equation (TlO|) does not allow wave splitting when evolutionary dynamics are decoupled 
from population dynamics. 



Effects of ecological on evolutionary dynamics 

To understand the effect of ecology on evolution, we assume that the ecological parameters do 
not depend on the frequency of defectors: K(f) = K, g c (f) = g c , and c*(f) = c*. The partial 
differential equation for c(t,x) is then independent from the dynamics of f(t,x), so the velocity 
of the population wave and wave profile are given by equations (J3J) and (jl]) respectively. As a 
result, equation (ITU]) is reduced to a single equation for the defector fraction f(t,x) by substituting 
the solution for c(t,x). In the reference frame (r, £) comoving with the population wave, this 
substitution results in an explicit dependence of the growth dynamics on ( through the c-dependence 
of Gf and 2£)JMn(c)|£. Compared to equation (J5J), the additional term 2Z}JMn(c)|^ comes from 
the effect of density gradients on spatial diffusion of alleles. This term is similar to advection in 
react ion- diffusion equations with the medium moving at a velocity v a = —2D-^\n(c). To make this 
analogy explicit, we need to take the advection velocity v a inside the partial derivative, which results 
in the following equation for f(r, C) 

% =D w + l [( " c ~ Va)f] + [Gf + ^ ]/ ' (12) 

where the term with v c appears because we changed to the comoving reference frame. From equa- 
tion f[T2"j) . one can see that 2D-^ ln(c)|^ results in both effective advection and effective growth. The 
effective terms are only functions of ( and can be computed from the known density profile c((). 
For c* > —K/2, this computation can be carried out explicitly with the following results 
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Figure 2: Splitting of cooperators from defectors. The top panels (a), (b), and (c) show a mixed 
wave where cooperators and defectors spread together. The bottom panels (d), (e), and (f) show 
cooperators splitting from defectors and colonizing new territories faster than defectors can invade 
them from behind. For all panels in this figure, we used K — 1, c* — 0.1/(1 — /), g c — 3, gf — 0.25, 
D = 0.5, and numerically solved equation (Q. The difference between the top and the bottom panels 
is in /*(c). No splitting occurs for /* = 5 $(c — 0.2) — 4.5, in the top panels, but a higher value of c 
ensures splitting in the bottom panels, where /* = 5$(c — 0.9) —4.5. (a) and (d) show a mixed wave 
expanding into new territories. In both panels, the front has a higher frequency of cooperators (blue) 
compared to the population bulk (magenta). However, the size of the region enriched in cooperators 
remains constant in the top panel, while it grows linearly with time in the bottom panel, (b) and (e) 
show the positions of the expansion fronts for cooperators and defectors as a function of time. The 
position of the front is defined as the rightmost point where the density of cooperators or defectors 
reaches half of its value in the bulk. The asymptotic rates of expansions are shown with dashed 
lines. From (e), we know the velocity of cooperators v c = 0.69 and the invasion velocity i\ = 0.44. 
These velocities are the same as in (b) because they do not depend on c. Note that the mixed 
velocity v m = 0.66 is smaller than v c because defectors increase the Allee threshold and slow down 
the population. More importantly, defectors, which are also spreading with velocity v m , colonize 
new territories faster than they can invade cooperators (v m > v{). (c) and (f) show density profiles 
of defectors and cooperators at t — 75. As expected, the lag between cooperators and defectors is 
larger in (f), where there is a region of pure cooperators at carrying capacity. 
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v a = -2D-?- ln(c) = ^lDg~K{\ - c/K) 

9e = "^T =9cC{K -C). 

The effective growth term is always positive and peaks at the middle of the population wave front 
where c = K/2 and ( = (by definition). The effective advection term is also positive, provided / 
decreases with £, which is expected when cooperators are favored at low population densities. As we 
show below and in the Appendix, these two terms in some cases allow defectors to keep pace with 
cooperators even when v\ < v c . 

The existence of a mixed traveling wave, where both cooperators and defectors spread at the 
same velocity, is equivalent to the existence of a steady state for f(r, £) in equation (fT2l) . When a 
steady state does not exist, population and genetic waves split, and f(r, () shifts to negative ( with 
velocity v c — f;. Since Gf decreases with /, the existence of a steady state requires that operator £, 
obtained by linearizing the right hand side of equation (fT2|) with respect to /, has a positive eigen- 
value A. Indeed, if all eigenvalues of C are negative, then ^ < because G/(c, /) < G/(c, 0), while 
a positive eigenvalue ensures that small /(r, () will increase until G/(c, /) is sufficiently diminished 
so that |£ = 0. Thus, we look for a solution of the following equation with A > 

Df"(Q + [v c - v.(C)]/'(C) + 21(0/(0 = A/(0, (14) 

where gi{() = Gf(c((), 0), and we used primes to denote derivatives with respect to (. Equation (|T4|) 
can be transformed to a canonical form by the following change of variables 

/(C) = e-^V(C) where 

which also insures that ip(±oo) = 0; see Appendix. The result reads 

- Dip" + 7(0^ = -A^>, where 

v(0 = - 9l (0 - i<(0 + ^K- ^a(C)] 2 - (16) 

We can now use the standard theory of second order differential equations to establish conditions 
necessary for the existence of a solution for A > 0; see [BH These conditions require that there 
must be a sufficiently large region where the potential V(() is negative and the values of V(() in 
this region must be sufficiently low. To show this, we multiply both sides of equation f lTB]) by ip and 
integrate over (, which gives 

/oo /»oo roo 

^(QdC + X ip 2 (()d(+ y(C)^ 2 (CR = (17) 
-oo J— oo J— oo 

after an integration by parts in the first term. Since the first two terms in equation ( ITT)) are positive, 
the third term must be negative, which in turn requires that there exists a region where V(Q < 
because i/j 2 (() > 0. This region has a finite width because V(±oo) > 0, which follows from 
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equations QS) and (III). Indeed, for C -> -oo, we find that V(-oo) = -g f (K)f*(K) + v 2 J(AD) = 
(v 2 — vf) / (AD) , which is positive since v c > v\. In agreement with the earlier discussion, when v c < v- u 
the potential at — oo is negative ensuring the existence of eigenfunction for A > and, therefore, 
the existence of mixed waves. For ( — > +00, we find that V(+oo) = — fl'/(0)/*(0) + (v c — f a ) 2 /(4D), 
which is positive as long as defectors are sufficiently disfavored at low densities, e.g. when /*(0) < 0. 
Thus, V(C) is a potential well, and the existence of an eigenfunction for A > requires that this 
potential well be sufficiently wide and deep. 

We now interpret the effects of the three terms on the right hand side of equation ( f!6l) in the 
context of the depth and width of the potential well V((). The first term lowers V(Q for c > c in 
the population bulk (£ — > —00), but it increases V(Q for c < c at the front (£ — > +00). The second 
term is always negative; it deepens the potential well around ( = and vanishes for £ — > ±00. 
The third term is always positive, but the reduction in the depth of the potential well due to v c is 
lessened by v a (at least for some (). The transition from non-splitting to splitting behavior can then 
be achieved by making the potential well more shallow. In particular, the transition threshold can 
be crossed by increasing c, the density at which natural selection starts to favor cooperators over 
defectors, as shown in figure |2j Decreasing /*(0) or v\/v c has a similar effect. 

To understand wave splitting better, we solve a special case of equation 014p exactly in the 
Appendix and find how the threshold between splitting and non-splitting behavior depends on model 
parameters. An exact solution can be found when c* = 0, and /*(c) = /* for c > c, but defectors 
are not viable for c < c. In other words, we impose an absorbing boundary condition at ( such 
that c(C) = c. With these simplifications, we find that splitting occurs provided 



We draw three important conclusions from this result. First, the exact solution confirms that both 
splitting and non-splitting behaviors are possible depending on the parameters of the model. Second, 
the severity of selection against defectors required for splitting increases with the ratio of invasion 
and cooperator velocities. Third, defectors can spread faster in mixed waves than they can invade 
cooperators when t> ; < v c , but splitting does not occur. These three conclusions do not depend on the 
simplifying assumptions used to derive equation ( fl8|) . which can be shown by numerical solutions of 
equation (jUJ) or ([16]) ; see, e.g., figure |2j We now discuss the biological significance of wave splitting. 

The possibility of wave splitting has important implications for the evolution of cooperation. 
When splitting is possible, cooperators outrun defectors, leading to a constant increase in the relative 
frequency of cooperators for as long as there are uncolonized territories. Frequent local extinctions 
followed by recolonizations could therefore maintain high levels of cooperation in natural popula- 
tions. Equation (JTSJ) also suggests that traits that have little effect on the dynamics in well-mixed 
populations might be under selection during range expansions. One such trait is gf. In well-mixed 
populations, it only determines the rate of approach to the equilibrium not the equilibrium itself, 
while, in expanding populations, gf affects the invasion velocity and the conditions for wave splitting, 
both of which are expected to be under selection. The critical density c is another examples. It has 
little effect in the well-mixed populations, but large c allows cooperators to escape from defectors 
during range expansions. 

The possibility of non-splitting highlights the importance of the coupling between ecology and 
evolution in the maintenance of genetic or species diversity during range expansions. Quite surpris- 
ingly, we find that defectors can spread in a mixed wave faster than they can invade cooperators 
despite the fact that defectors are disfavored or even eliminated by natural selection at the front. 
In other words, the community of cooperators and defectors is able to move to new territories while 
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preserving the coexistence between the two phenotypes. Splitting, on the other hand, would lead to 
a loss of the community structure. Our results could, therefore, have important implications for con- 
servation efforts to preserve genetic diversity or ecosystem integrity during potentially rapid range 
shifts due to climate change or habitat deterioration. For negative frequency-dependence discussed 
here, we found that diversity is more stable than one would naively predict from just measuring v c 
and v\. Other types of interactions could be less resilient and would require managed interventions 
to prevent splitting. In the Appendix, we show that interventions increasing advection or growth 
rate at the front can achieve that goal. 

General case of eco-evolutionary feedback 

In the general case when ecology and evolution affect each other, the nature of the transition from 
non-splitting to splitting behavior remains the same. In particular, the condition for splitting is 
still given by the existence of a solution of equation ( Fl4l) with A > because, close to the splitting 
transition, the expansion front is populated almost exclusively by cooperators. To show this, let us 
consider how the dynamics at the front changes as one of the model parameters, say c, is varied 
so that the system behavior changes from non-splitting to splitting. Increasing c directly increases 
the abundance of cooperators at the front and increases the velocity of the mixed wave from Vf*, 
when c = 0, to v c , when splitting occurs. In addition to that, defectors lag more and more behind the 
front, as the splitting transition is approached. Indeed, close to the transition, Gf is barely sufficient 
for defectors to keep up with cooperators, and, since Gf decreases with /, the fraction of defectors 
at the front must approach zero right before splitting occurs. These dynamics are illustrated in 
figure [3^ and further discussed in the Appendix. Another effect of the coupling between ecology and 
evolution is that V\ is no longer given by equation (J7|) and one has to solve equation (fit)]) to describe 
defector invasion. 

We would now like to discuss the effect of initial conditions on wave splitting. One may naively 
think that deterministic wave splitting can be achieved simply by creating a region of pure coopera- 
tors in front of a mixed population, even if populations do not split when started from a well-mixed 
state. Initially, the range expansion started from such an initial condition indeed resembles splitting 
with cooperators transiently outrunning defectors and increasing in numbers; see figure [3b. Never- 
theless, defectors eventually catch up because any initial condition (with c(0, x) > and f(0,x) > 0) 
has a nonzero projection on the eigenfunction of the operator C with the largest positive eigenvalue. 
The growth of this projection ensures the establishment of defector population at the front even 
when v\ < v c as shown in figure [3b. The ability of defectors to catch up, however, relies on the 
assumption that c(t,x) and f(t,x) vary continuously and can increase from arbitrarily small values. 
When the discrete nature of organisms is taken into account, one finds that wave fronts have a finite 
width [HS]- Therefore, initial conditions can force splitting provided v c > v- x and the region of pure 
cooperators is sufficiently large compared to the width of expansion and invasion wave profiles. 

Number fluctuations must also be considered for traveling waves of discrete entities. For bi- 
ological systems, these fluctuations could arise due to random fluctuations of the environment or 
due to the randomness of births and deaths, which is known as genetic or ecological drift. Number 
fluctuations are largely irrelevant when V[ > v c because defectors can always catch up with cooper- 
ators. In contrast, when v\ < v c , splitting is the ultimate outcome because, given enough time, a 
fluctuation will create a region of pure cooperators, which is sufficiently large to prevent defectors 
from catching up. This region of pure cooperators will then grow with time making it exceedingly 
unlikely for another fluctuation to destroy it. Although splitting is inevitable, it may take a very long 
time (e.g. 10 4 generations in figure [3b) because the deterministic dynamics discussed above create 
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Figure 3: Wave splitting under different conditions, (a) shows how the velocities of cooperators (blue 
circles), defectors (red dots) depend on c. For small, c the population moves with the velocity Vf*. 
As c increases, the expansion first reaches velocity v c and then splits into a population and a genetic 
wave, which occurs when the velocities of cooperators and defectors are no longer the same. The 
velocity of cooperators and the population expansion is v c , while the velocity of defectors is V{. For 
this panel, we used K = l,g c = 2,g f = 0.05, c*(/) = 0.24/(1-/), /* = lO0(c-c)-9.5, and D = 0.5. 
We used large habitat lengths up to 1800 because, close to the splitting transition, transient dynamics 
decay slowly, (b) Creation of a region of pure cooperators in front of a mixed population leads to 
transient splitting. Initially, cooperators seem to outrun defectors because v c > v\. At later times, 
however, defectors catch up because the conditions for deterministic splitting are not met, and a 
small number of defectors that migrated to the front can follow cooperators and grow until a stable 
mixed wave is formed. For this panel, we used K — 1, g c — 3, gf — 0.25, c*(/) = 0.1/(1 — /), 
/* = 5$(c — 0.1) — 4.5, and D = 0.5. (c) Stochastic splitting is possible even when the conditions 
for deterministic splitting are not satisfied, but v c > V{. The waiting time for stochastic splitting can 
however be very long (here 10 4 generations), even for moderately low population densities of 100 
individuals per site used in this simulation, because splitting is caused by very rare fluctuations that 
create a large region of pure cooperators at the front. For this panel we used L = 1000, m = 0.1, 
K = 100, g c = 0.05, c* = 0, and the entries of the interaction matrix are given by A cc = 2.5, 
A c & = 1.5, Adc = 2.65 and = 1. 



an effective activation barrier for stochastic splitting; see figure [3b and HI The magnitude of this 
barrier goes to zero as the conditions for deterministic splitting are approached. Stochastic splitting 
is also possible for frequency-independent selection, provided expansion and invasion velocities are 
different; see for example |42j . 

In the preceding discussion, we neglected possible transitions between cooperators and defectors 
due to mutations or other heritable changes. This is justified on short time scales because mutation 
rates are typically small and even beneficial mutations struggle to survive number fluctuations (ge- 
netic drift) [67J ES]. On long time scales, mutations will change the dynamics in two ways. First, 
the coexistence fraction /* will be determined by both natural selection and the mutational pres- 
sures similar to the classic theory of two genetic variants [6TI 133] . This shift in /* does not change 
the dynamics qualitatively and could be included in our theory by modifying Gf(c, /) accordingly. 
The second and more important difference is that the region of pure cooperators will not expand 
indefinitely following a splitting event because defectors will appear in the interior of the region of 
pure cooperation due to a mutation in addition to invading the region of pure cooperators from 
behind. In the limit of rare mutations, the average length of the region of pure cooperators L pc will 
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Figure 4: The average time to stochastic splitting shows an exponential-like increase with the car- 
rying capacity of the sites K. The parameters used in these simulations are the same as in figure [3fc 
with the exception of L = 2000 for K < 100 and L = 3000 for K = 100. The error bars show the 
standard deviation of the mean. Similar increase in the average time to splitting was observed for 
increasing migration because larger m leads to longer wave profiles containing more organisms, and, 
therefore, smaller number fluctuations. 



be determined by the balance between the time necessary to create this region after splitting and 
the time to the next successful defector mutation in the region of pure cooperators. The former time 
scales as L pc /(v c — v{) for deterministic splitting. While the latter time scales as the inverse of the 
product of the mutation rate from cooperators to defectors fid, population size K(0)L pc , and fixation 
probability g f [K(0)]f*[K(0)], i.e. as {/id Lp^^)^^^)]/*^^)]}- 1 ; see [6T 1[35]. For determinis- 
tic splitting, this balance yields L pc m (v c — Vi)/{fj,dK(0)gf[K(0)]f*[K(0)]}, while, for stochastic 
splitting, L pc will also be affected by the average waiting time to stochastic splitting. Note that, even 
when a mutation appears right at the expansion front, it would not lead to a permanent collapse 
of the region of pure cooperators when mixed defector-cooperator waves split deterministically. In 
the other limit of very frequent mutations, splitting would not be able to occur before additional 
defectors arise at the front, and the dynamics would be primarily determined by the balance between 
mutational transitions leading to mixed populations. 

Another interesting consequence of the coupling between ecological and evolutionary dynamics 
is wave acceleration. We found that cooperators are favored at the wave front (see figured]). As 
the frequency of cooperators is increasing at the front, the instantaneous velocity of the expansion 
wave must also increase because defectors slow down expansions (see figure [2]). The evolutionary 
change could however be very slow compared to the colonization dynamics leading to a gradual 
acceleration of the range expansion. Indeed, long periods of dramatic wave acceleration are possible 
in our model and are shown in figure The acceleration of expansion fronts has been observed in 
a number of species and is typically explained by the evolution of new traits like shorter generation 
times, greater dispersal abilities, or better adaptation to the environment in the newly colonized 
regions [681 |3T1 [691 [30] . Our analysis of cooperator and defector waves suggests that a changing rate 
of expansion does not require the evolution of new traits, but could simply be a consequence of the 
slow adjustment of the genetic composition at the wave front to a new low-density optimum, which 
is different from that in the population bulk. 
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Figure 5: Population waves can accelerate due to evolutionary changes at the front, (a) A mixed 
wave of cooperators and defectors spreads into empty territories. The shape of this invasion is not 
linear in the (t, x) coordinates indicating wave acceleration. The acceleration is caused by a slow 
increase in the number of cooperators (blue) at the front compared to the bulk (magenta), (b) 
The velocity of the front increases from Vf* to its asymptotic value v m , which is smaller than v c 
because population and genetic waves do not split in this case, (c) Cooperator and defector density 
profiles at different times show both acceleration of the expansion and the increase in the frequency 
of cooperators at the front. For this figure, we used K — 1, g c — 5, gf — 0.1, c*(f) = 0.24/(1 — /), 
/* = 3 (c - 0.5) - 2.5, and D = 0.5. 



Discussion 

Understanding the link between ecology and evolution remains a major open problem [70]. This 
coupling is particularly important during range expansions, which are often accompanied by both 
demographic and evolutionary changes [2H EHJ EH EE]- We formulated a simple two-allele model 
that is capable of describing a wide range of frequency and density dependencies of selection and 
growth. Our analysis revealed the necessity of taking eco-evolutionary feedback into account in order 
to make accurate predictions. In particular, we found that the genetic composition of the population 
bulk may be a poor predictor for the genetic composition of the population front. Similarly, the 
measurements at the expansion front may be a poor predictor of the properties of the new population 
once it is fully established. These differences between the population bulk and expansion front make 
it also harder to predict the rates of expansions. Indeed, we found that the rates of expansions could 
be accelerating as the result of slow evolutionary changes at the expanding population edge and 
that species can spread faster in mixed waves than in isolation or by invading already established 
populations. 

Our main result is that colonization of new territories can proceed in two qualitatively different 
ways. The first possibility is a single mixed wave, where both alleles (or species) move with the same 
velocity and their relative abundances reach a steady state in the reference frame comoving with the 
expansion. The second possibility is that one of the alleles (or species) outruns the other. The faster 
allele is solely responsible for the colonization, while the slower allele invades the faster one from 
behind with a smaller velocity. As a result, there is a growing region occupied exclusively by the 
faster allele. The effect of wave splitting could be especially dramatic for species that have markedly 
smaller migration rates in the population bulk compared to population front [7TJ H2] because the 
invasion by the slower allele would be significantly slower than the range expansion. The existence of 
secondary genetic waves could lead to unexpected population and genetic dynamics, which cannot 
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be described by the classic models of range expansions [48l 1471 [571 [56] . It would be interesting 
to know when the commonly observed changes behind the expansion front [69] are caused by de 
novo adaptation to the new environment and when they are caused by the secondary genetic waves, 
which could reach and alter populations at the newly colonized territories many generations after 
the arrival of the population wave. 

In the context of cooperator-defector interactions, wave splitting and the increase in cooperation 
at the front could stabilize cooperation against defectors in populations that experience frequent 
disturbances followed by recolonizations. This mechanism of maintaining cooperation does not rely 
on reciprocity or multi-level selection, but is instead grounded in the density dependence of the 
evolutionary dynamics. More importantly, we have demonstrated that the coupling between ecolog- 
ical and evolutionary dynamics can have a profound effect on the fate of cooperation and should, 
therefore, be considered in both theoretical and experimental studies. 
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Appendix 

In this Appendix, we first expand our discussion of the boundary conditions for equation (fT6|) . We 
then outline the derivation of the exact solution, given by equation ( JT8l) . for the transition from 
non-splitting to splitting behavior in a special case of equation ( fT2l) . Finally, we show that both 
effective advection and effective growth can allow defectors to expand faster in mixed waves than 
they can do by invading cooperators. 



Boundary conditions for equation ( T76l) 

Here, we show that ^(±oo) = are the proper boundary conditions for equation ( jTBT) . Since the 
boundary conditions for f(Q are /(— oo) = 1 and /(+oo) = 0, it is sufficient to show that u(— oo) > 
and u(+oo) < 0; see equation (fT5|) . The former condition follows from the fact that v c > 
and v a (— oo) = 0. The latter condition follows from the fact that v a (+oo) > v c , which we prove below. 
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The density profile at large ( decays exponentially as e~ 7< >; moreover, 7 > ^/G c (0, 0)/D for a 
weak Allee effect [73]. We then substitute this scaling form in the equation ([9]) for c(t, x), change 
to the comoving reference frame, and require that the wave profile does not depend on r. This 
procedure yields 

v ^ Dl + ^Ml. (19) 

7 

This expression for v c can now be compared to the advection velocity v & = — 2Ddln[c(t, x)]/dx = 
2D7. For a strong Allee effect, G c (0, 0) is negative and, therefore, v a > v c . For a weak Allee ef- 
fect, v a — v c = D^f — Gc ^ ) ' - > > 0, where the last inequality follows from 7 > ^G c (0, 0)/D; see 



Derivation of equation / flgj) 

Here, we derive equation ffTgj) . Since, for Gf(c,f) decreasing with /, wave splitting is affected 
only by the linearization of Gf(c,f)f for small /, we are free to replace Gf(c,f)f with any other 
function that has identical behavior for / — > 0. It is particularly useful to approximate Gf(c,f)f 
with a piece-wise linear function because this allows for an exact solution for f(t, x) everywhere in 
space, which captures the qualitative dependence of the solution on the parameters of the model. In 
contrast, eigenf unctions of £ only describe the behavior of f(t,x) close to the front and only up to 
a normalization factor. To this purpose, we redefine G(c, f)f as 

G f (c, /)/ = g f (c)f*(c) min{/, f (c) - /}. (20) 

Note that this redefinition does not have a fixed point at / = 1, which is appropriate as long as we 
are not interested in cooperators invading defectors. 

Equation ( TT2|) can be further simplified by the change of variables from ( to tp = c/K using 
equation (j4j). The result reads 

^(1 _ tff + 2^(1 - <p) (l - 2<p + £\ + -^[G f (Kp, f) - A]/ = 0, (21) 

where primes now denote derivatives with respect to tp. We further assume that c* = and obtain 
that 

b 2 (l - V) 2 H = --^G f (K<p, f)f, (22) 

where we also set A = because A vanishes at the splitting transition. We now use equation (j20l) 
and the assumption that f*(<p) = f* for ip > (p = c/K and f(tp) = for p < <p, i.e. defectors are not 
viable at very low densities. The right hand side of equation (1221) is then proportional to / for / < 
/*/2 and to /* — / for / > /*/2; therefore, further simplifications occur upon integrating both sides 
of equation fl22l) over ip and defining a new variable F — J fdtp for / < f*/2 and E = J (/* — f)dtp 
for / > /*/2. For example, for the region where / < f*/2 and <p > tp, we obtain 

v\l-v) 2 F" = -m^F (23) 
The general solution of this equation is given by 
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F = Cnp a {l - ip) 1 - + CV~ Q (1 - p) a , with 

8 9 ff* \ (24) 




Similarly, for / > /*/2, we find 

E = C 3 /(l - ^) 1_/3 + C 4 ^ 1-/3 (l - with 




(25) 



Then, the solution for f((p) is given by E' from equation (1251) for / > /*/2 and by F' from 
equation ( l24|) for / < f*/2 and tp > <p. The four constants Ci, C2, C3, and C4 and the matching 
point (^1/2 are determined by the boundary conditions /(— 00) = /* and f(c/K) = and the 
matching conditions at <pi/2, where /(</?i/2 — 0) = /(</?i/2 + 0) = 1/2 and f'(<fi/2 — 0) = f'{<pi/2 + 0). 
Upon using four of the five conditions, one can easily express the four constants in terms of if 1/2 
because the conditions are linear with respect to C±, C2, C3, and C4. The remaining condition 
determines y?i/ 2 and can also be solved exactly: 



V1/2 



7 



7^ 



where 



l-(l-7)^ 
(g - q)(q - g) 
_(/? + « -1)(1 -a -£) 



(26) 



We then obtain the condition for splitting, stated in equation ( {TBI) , by requiring that 7 is real, 
and ifx/2 > (p. Note that, ipi/2 -> 1 as the splitting state is approached, say, by increasing c. 
Since ip — > 1 corresponds to £ — >■ —00, we see that the wave of defectors lags more and more behind 
the wave of cooperators. As we show in the main text, this is a general result, and the front is pop- 
ulated mostly by the cooperators close to the splitting transition. As a consequence, the expansion 
velocity at the splitting transition is given by equation flS}. 



Increased spreading through effective advection and growth 

Here, we demonstrate the effects of effective advection and growth on the ability of defectors to keep 
up with cooperators. When ecological dynamics are not affected by the evolutionary dynamics, the 
ability of defectors to follow cooperators in a mixed wave is mathematically equivalent to the ability 
of a Fisher wave to follow a moving boundary between good and bad conditions. In particular, 
the assumption that / = for c < c makes this moving boundary absorbing. One can easily show 
that classic Fisher waves can only follow an absorbing boundary that moves with a velocity smaller 
than the Fisher velocity. Here, we show that even boundaries moving with higher velocities can 
be followed provided there is an effective advection or growth. This analysis illustrates the role of 
the term 2D 9 1 ^ |^ in equation (1101) and suggests that both effective advection and growth could 
be useful management strategies to maintain genetic diversity during range expansions or to help 
species shift in response to a rapid climate change. 

To avoid confusion with the earlier discussion, we introduce new notation. Let v be the velocity 
of the absorbing barrier, and g is the per capita growth rate. The Fisher velocity is then 2y/gD, 
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where D is the effective diffusion constant. When v > 2y/gD, the population cannot keep up with 
the barrier without external interventions. One possible intervention is to increase g by a factor 1+8, 
say by supplying extra food, in a region of length a behind the barrier, which is analogous to the 
effective growth discussed earlier. Another possibility, analogous to effective advection, is to move the 
population with velocity u in a region of length a behind the front, say by capturing and transporting 
individuals within the population. The analysis of both of these scenarios follows the same procedure 
as the derivation of ( IT8l) . but is much simpler because the resulting equations are linear and have 
constant coefficients; therefore they can be easily solved by the standard methods. Below we just 
summarize the results. 

For effective growth, we find that barriers moving with any velocity can be followed, but the 
region of effective growth should be sufficiently large: 



a > 



2D 



y/4gD(l + 6) - v 2 



71 



arctan 



/ y/4gD(l + 6)-v* 
\ v^ 2 - ^gD 



(27) 



Note that the minimal size is finite when v = 2y/gD, and that it decreases only as <5 _1//2 for large S. 
This suggests that there is an optimal size a that minimizes 8a, which could be interpreted as the 
total cost of the intervention. 

For effective advection, we also find that barriers moving with any velocity can be followed 
provided 2^gD + u > v and the region of advection is sufficiently large: 



a > 



2D 



^AgD -{v- uf 



tx — arctan 



y/4gD-{v-^ 
u + Jv 2 - AgD 
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Similar to the previous case, minimal a is finite when v = 2y/gD. More importantly, the minimal 
size is the smallest when u and v are about equal. Indeed, when v ^> u, advection has a small 
effect on the dynamics, while, when u ^> v, the point where the region of advection ends becomes 
an effectively absorbing boundary because organisms entering advection region are quickly moved 
towards the front of the wave. 
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